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Abstract 

We study the effect of gradual symmetry breaking in a non-integrable 
system on the level fluctuation statistics. We consider the case when 
the symmetry is represented by a quantum number that takes one of two 
possible values, so that the unperturbed system has a spectrum composed 
of two independent sequences. When symmetry-breaking perturbation is 
represented by a random matrix with an adjustable strength, the shape 
of the spectrum monotonously evolves towards the Wigner distribution as 
the strength parameter increases. This contradicts the observed behaviour 
of the acoustic resonance spectra in quartz blocks during the breaking of a 
point-group symmetry that has two eigenvalues, where the system changes 
in the beginning towards the Poisson statistics then turns back to the 
GOE statistics. This behaviour is explained by assuming that the sym- 
metry breaking perturbation removes the degeneracy of a limited number 
of levels, thus creating a third chaotic sequence. As symmetry breaking 
increases, the new sequence grows at the expense of the initial pair until 
it overwhelms the whole spectrum when the symmetry completely disap- 
pears. The calculated spacing distribution and spectral rigidity are able 
to describe the evolution of the observed acoustic resonance spectra. 

PACS: 05.45. -Kb, 11.30. Er, 24.60.Lz, 62.30+d 
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1 Introduction 



The statistical theory of spectra [Q provides an appropriate method for exam- 
ining the symmetry properties of quantum systems. Level statistics such as the 
nearest-neighbour spacing (NNS) distribution and the spectral rigidity A3 are 
defined for a pure sequence of levels that have same quantum numbers. Usually, 
it is taken as granted that in a classically integrable system, which has as many 
integrals of motion (or quantum numbers) as the number of degrees of freedom, 
the levels are uncorrelated and so have a Poissonian NNS distribution. The 
pure level sequence in a time-reversal-invariant quantum system whose classi- 
cal counterpart is chaotic, is successfully represented by a Gaussian orthogonal 
ensemble (GOE) of random matrices. A deviation of level statistics from the 
predictions of the GOE may be taken as an evidence that the sequence is not 
pure because of the existence of a hidden symmetry. The statistical theory of 
spectra has also been extensively used in the investigation of symmetry break- 
ing because the destruction of a quantum number has a dramatic effect on the 
spectral fluctuations. Impressive are the studies of isospin mixing in light nuclei 
Shriner et al.^ measured a "complete" spectrum of low- lying states of ^^Al 
involving states with isospins T — and 1 throughout the energy range covered 
by the data. These data offered a testing ground for studying the influence 
of isospin-symmetry breaking on the fluctuation properties of energy spectra 
(see, e.g., |4[ ||, ||, ^). However, the limited data available in ^^Al (142 levels) 
precluded a definite conclusion. In a recent experiment with monocrystalline 
quartz blocks, EUegaard et al. measured about 1400 well-resolved acoustic 
resonances. The properties of quartz allowed these authors to measure the grad- 
ual break-up of a point-group symmetry, which is statistically fully equivalent 
to the breaking of a quantum number like isospin. The spectra are obtained by 
externally tuning the symmetry breaking Q , allowing a statistically significant 
investigation of the whole transition. 

In the present paper we consider the effect of gradual symmetry breaking 
on the fluctuation properties of energy levels. When the symmetry is full the 
spectrum is divided into independent sequences, each corresponding to one of the 
eigenvalues of the symmetry operator. We consider two possible scenarios for 
breaking the symmetry. The first is provided by the Guhr-Weidenmiiller model 
, which assumes that the perturbation acts equally on all states and that the 
matrix elements between states of different symmetry eigenvalues are random 
numbers of equal variances. In the second scenario, the perturbation acts 
differently on the unperturbed states and the symmetry breaking interaction 
mixes a limited number of the degenerate eigenstates that belong to different 
symmetry representations, and this number grows as we increase symmetry 
breaking. We are thus creating a new sequence of levels corresponding to states 
with no symmetry at the expense of the initial sequences. Consequently, as the 
fractional density of the no-symmetry sequence increases, the A3 statistics of 
the total spectrum first increases, moving in the direction of the prediction of 
the Poisson statistics until the densities of all the sequences become equal, and 
then decreases to reach asymptotically the value given by the GOE. Section 3 
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shows that the behaviour described by the second scenario is indeed present in 
the data of EUegaard et al. Q . The summary and conclusions of this work are 
given in Sect. 4. 



2 Symmetry breaking in a chaotic system 



We consider a chaotic system that conserves a given symmetry. In such a 
system, the Hamihonian assumes a block-diagonal form 



i/5 = diag {Hi,H2,---,H.n), 



(1) 



where each sub-block corresponds to states of certain quantum numbers. Thus, 
when the symmetry is conserved, the levels are divided into a number of "pure" 
sequences, each described by a pure GOE. However, the statistical properties 
of the total spectrum that combines these sequences are no longer of the GOE 
type. In the limit of infinite number of sequences combined, the spectrum 
obeys the Poisson statistics. For a finite number of sequences, the fluctuation 
properties are intermediate between the Poisson and GOE statistics. In the fol- 
lowing, we shall restrict our consideration to the case when the symmetry under 
investigation is represented by an operator that has two possible eigenvalues, 
such as parity, or isospin and 1 as considered in Refs ^, H, ||, Q . 

The NNS distribution of a spectrum resulting from a random superposition 
of n independent sequences is calculated, e.g., by Berry and Robnik ||^ and in 
Mehta's book . If the level density of the ith sequence is pi , and if the NNS 
distribution of levels of this sequence is Pi{xi), where Xi — fis, fi — pifSpi, 
and s is the NNS normalised to a unit mean, then the NNS distribution of the 
mixed sequence is given by 



H 2 



E 



(2) 



where E{s) = 11^=1 E^[Us), E,{x,) = [1 - (x)] dx, and W, (xi) = J^' P,{x)dx. 
In particular, if all the n individual sequences have the same level densities, so 
that fi = 1/n, and if the NNS distribution in each is a Wigner distribution 



Piix^) = -x^e i'- 



(3) 



which is a good approximation for the NNS distribution of a GOE, then (2) 
becomes llli 



Puis) 



erfc 



2n 



(s) [^s + (n-l)Q„(s) 



(4) 



where Qn{s) — exp (— 7r^s^/4n^) /eik{sy/TT/2n) and erfc(a:) is the complemen- 
tary error function. This distribution has a shape intermediate between those 
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of the Wigner and the Poisson distributions. The Brody distribution 

Pe (.) = e,.«exp (-^.-') with c, ^ E!!lM£±i)), (5) 

accurately reproduces the spectra of many two-dimensional systems with mixed 
regular-chaotic dynamics. As we vary the parameter /3 from to 1 the statistics 
changes from Poissonian to GOE. 

The NNS distribution contains information about the spectrum in a short 
range, not exceeding a maximum of three mean level spacings. Long range 
information is provided by various higher-order correlation functions p^ . 
The most popular among these is the spectral rigidity A3 of Dyson and Mehta 
[ p!3| . Semiclassical arguments |Q show that, for a generic integrable system 
whose spectrum satisfies Poisson statistics 

A3,Poisson(^) = ^ for L < Lmax, (6) 

where imax = 27r7i/DTniin, D is the mean level spacing and T,„iii is the period 
of the shortest classical orbit. For time-reversal-invariant chaotic system the 
semiclassical theory gives the following asymptotic expression 

A3,GOE(i) = ^ In L- 0.007 for 1 « L « (7) 

In all cases A3 saturates to a non- universal value at L ^ ^max- The random 
matrix theory obtains an analytical expression for A3_GOE(i) that involves dou- 
ble integration |l^ . We find it more suitable to parametrise this expression 
for a GOE, guided by Eq. (0), in the form 

A3.GOE(i) - ^ (1 - e""'') ML + b) + c], (8) 

with a — 8.241,6 — 0.944 and c = —0.06724 being the best fit parameters for 
values of L in the range of < L < 50. Seligman and Verbaarschot sug- 
gested the following expression for the spectral rigidity of a spectrum resulting 
from a random superposition of n independent sequences with fractional level 
densities fi 

n 

A3(L) = ^[A3,.(/.i)], (9) 

1=1 

where A3.i(L) is the spectral rigidity for the ith sequence. This proposal is not 
surprising because A3(L) is essentially a variance. 

We now consider symmetry breaking in the chaotic system using the Guhr- 
Weidenmiiller model The model suggests that the effect of the symmetry- 
breaking perturbation on all the unperturbed states is equal on the average. 
This suggests defining the Hamiltonian in the form: 



H 



' Ho 










w ' 





Ho . 


+ a 








(10) 
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where Hq represents a GOE and Vl^ is a random matrix whose elements have 
Gaussian distributions with zero means and variances equal to those of the non- 
diagonal elements of Hq. When a = 0, the symmetry is conserved and the 
Hamiltonian has the block diagonal form of Eq. (|^) with n — 2. The spectrum 
consists of an independent superposition of two level sequences having equal 
density. The NNS distribution P(s) is given by Eq. and, in particular, 
P(0) = 1/2. The spectral rigidity is given by (|) with /i = /a = 1/2. When 
we allow a to increase, the system gradually evolves toward a single GOE. To 
show how this transition proceeds, we considered an ensemble of ten 200x200 
Hamiltonian matrices of the form ( |l0[ ) for a fixed value of the perturbation 
strength a. We numerically diagonalised these matrices. We obtained the 
NNS distribution and spectral rigidity for each of the resulting spectra, and 
measured their deviation from the GOE statistics by comparing with the Brody 
distribution, Eq. (|^), for the NNS distribution and with the Seligman and 
Verbaarschot expression, Eq. (||), for the A3 statistics where Aa^i is given by 
the Poisson distribution (g) and A3^2 by the GOE (ph, so that 



A, 



3,SV 



qsv— + A3, GOE [(1 - qsY)L] , 

15 



(11) 



We then repeated the procedure for other value of a. The best fit of the tuning 
parameters f3 and qsv obtained for P{s) and A3(i) corresponding to different 
values of the perturbation strength a are shown in Fig. 1. The error bars 
measure the changes of the parameters that increase the minimum values 
by a factor of 2. The monotonous nature of the curves clearly shows that the 
transition from symmetry to no-symmetry in the Guhr-Weidenmiiller model is 
gradual and one-directional. 

The opposite extreme of the Guhr-Weidenmiiller model is the case when the 
eigenfunctions will become either very weakly or very strongly perturbed. For 
those states in which the perturbation is weak enough, the wave function can 
in principle be characterised by one of the two symmetry eigenvalues. These 
will be the states for which the splitting of two degenerate levels belonging to 
the two symmetry representations can be neglected. The corresponding energy 
levels will again form two independent sequences with fractional level densities 
/i = /2 = //2, while the levels of the strongly perturbed eigenstates will form a 
new sequence with density — 1 — / which will grow at the expense of the other 
two as the symmetry-breaking perturbation increases. The NNS distribution 
will then be given, using (||) and @, by 



P(5) = |/3,e-/^^Vi6e,fc(^\^/s^erfc 
+ i/^e--^^^V«erfc[^(l-/). 



+ |(l_/)3,e-.(i-/)^sV4 



erfc — — fs 
4 
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+2/(1 - /)e--(5/'-«/+4)^Vi6e,fc (^^/s) . (12) 

for spacings greater than a certain lower limit. If we ignore the latter condition, 
the extrapolation of P{s) into the origin yields 

F(0) = /(^2-^/), (13) 

We note that P(0) increases with decreasing / from a value of 1/2 at / = 1, 
reaches a maximum value of 2/3 at / = 2/3, and then decreases monotonously 
until it vanishes at / = 0. Equations (12) and ( [l^ ) show that the evolution 
of the shape of the NNS distributions with increasing symmetry breaking, and 
decreasing /, is not a straightforward transition from the 2-GOE behaviour to 
that of a GOE. Indeed, by decreasing the value of / we are removing some of 
the levels from the initial two symmetry-invariant sequences to form the third 
sequence that represents the states without symmetry. The distribution will 
look more regular when we decrease / until the level densities of the three 
sequences become equal. If we further increase the symmetry breaking, the 
initial sequences become thinner until the whole spectrum forms a single GOE 
sequence. 

We can follow this transition better by calculating the spectral rigidity which 
is, for this system, given by 

A3 (L) = 2A3,GOE (^^) + As.GOE [(1 -f)L]. (14) 

In Fig. 2, the spectral rigidities corresponding to / = 1, 0.9, 0.667, 0.333, 0.1 
and are given by the curves labelled by A, B, C, D, E and F, respectively, 
while the curve P corresponds to a regular spectrum. The figure clearly shows 
that A3 (L) increases as symmetry breaking increases, and thus evolves towards 
the regular shape until / reaches the value of 2/3. Only then starts A3 (L) a 
route leading to the GOE curve. 



3 Comparison with acoustic resonance spectra 

In this section, we show that the proposed application of the Berry-Robnik 
theory to symmetry breaking leads to a satisfactory description of the acoustic 
resonance spectra of monocrystalline quartz blocks, measured by EUegaard et al. 
[||. Crystalline quartz exhibits D3 point-group symmetry about the crystal's 
Z (optical) axis, and three two- fold rotation symmetries about the three X 
(piezoelectric) axes; the latter three axes lie in a plane orthogonal to the Z axis 
and sustain angles of 120 degrees with respect to each other. EUegaard et al. 
used rectangular blocks of dimensions 14 x 25 x 40 mm^ , cut in such a way that 
all symmetries are fully broken except a two- fold "flip" symmetry about one 
of the three X axes. A gradual breaking of this symmetry is then achieved 
by removing successively bigger octants of a sphere from one corner, thereby 
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creating a three-dimensional Sinai billiard. The acoustic spectrum is measured 
for 8 radii of the removed octant, providing data of high statistical significance 
for the study of the gradual breaking of a point-group symmetry. 

Figure 3 shows the spectral rigidity for the cases when the symmetry is 
present (open circles), and when it is violated by removing a tiny octant of 
radius r = 0.5 mm (closed circles). The straight line labelled P is for an inte- 
grable system. We note that data for the case of tiny symmetry violation are 
systematically higher than the ones for the case of conserved symmetry. This 
agrees with the prediction of the second scenario proposed above, in which the 
system that undergoes a gradual symmetry breaking looks as if it has become 
more regular in the early stages of the transition. 

EUegaard et al. concluded from the rise of A3 over the theoretical pre- 
diction for a superposition of two GOE sequences that the quartz block with 
conserved flip symmetry has much in common with a pseudo-integrable sys- 
tem. They supported this conclusion by means of other measurement. Pseudo- 
integrable (PI) systems are non-integrable systems, yet non-chaotic. Two- 
dimensional PI systems, exemplified by particles moving in a planar polygonal 
enclosure with rational angles [^, share with the integrable system that the 
motion of trajectories in the phase space is restricted to two-dimensional com- 
pact surfaces. But these invariant surfaces are multi-connected and not tori as 
in the case of integrable systems. Numerical studies of the 7r/3 rhombus bil- 
liard [|l^, |l^, |l^ show that the spectral fiuctuation properties are intermediate 
between those of a Poisson ensemble and a GOE. Bogomolny et al. Q studied 
the level statistics of PI systems using a modified version of Dyson's stochastic 
Coulomb gas model , in which the interaction between particles is restricted 
only to the nearest neighbours. The NNS distribution obtained in is given 

by 

Ppi(s) =4se-2s^ (15) 

This distribution is in agreement with the numerical calculation for the rational 
billiards ||l8|. Unfortunately, the formula obtained by these authors for the 
two-point correlation function 

i?2(s) = l-4e-4« (16) 

leads to an expression for the A3 statistics that asymptotically increases as 
L/30, which is not consistent with the results of the numerical experiments. 
Following several authors, e.g. |lj], we shall represent the A3 statistics by 
Seligman-Verbaarschot ||l^ formula (||), obtained for a superposition of Poisson 
and GOE level sequences, 

A3.Pi(L) - -7-^ + A3.G0E [(1 - q)L] , (17) 
15 

where q is a fitting parameter. Biswas and Jain found q to be equal to 0.2 
for the 7r/3 rhombus billiard. 

We now apply the second scenario, which has been proposed in the previous 
section, to the symmetry breaking in a PI system having a symmetry repre- 
sented by a quantum number that takes one of two possible values. When 
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the symmetry is conserved, the spectrum is an independent superposition of 
two sequences corresponding to the two representations of the symmetry. The 
NNS distribution is accordingly represented by (||) and the A3 statistics by 
with n — 2, Pi (xi) given by ( p^ ) and A3_i(L) by (p7|). We now again assume 
that as the symmetry violation starts, the eigenfunctions of the system will be 
divided into two classes. The first constitutes a fraction / of the eigenfunc- 
tions in which the perturbation due to symmetry breaking can be ignored. The 
levels of these eigenstates are divided into two independent sequences exactly 
as before switching the symmetry breaking on. The second is the class of the 
strongly perturbed states that have completely lost the symmetry. This latter 
class will here be modelled by a GOE. Therefore, the total spectrum will be 
composed of 3 independent sequences, two with PI statistics and each having 
a fractional density of //2, and one with GOE statistics and fractional density 
1 — /. During the symmetry-breaking transition, the NNS distribution is given 

by 



P{s) = (2/2^2 + 4/s + l)e-2/-erfc 



+ - /)(2 + fs) X [iTfs'il - ff + 2s {(tt + 4) f - 2nf + tt} + 8/] 



which is obtained by substituting Eqs. (^) and (15) into Eq. The A3 is 

given by 

A3 {L) ^ 2A3.P1 (^^^ + A3,GOE [(1 -f)L]. (19) 

We now compare our model with the acoustic resonance spectra measured by 
Ellegaard et al. We start by estimating the parameter q in the expression 
( p7| ) for the A3 statistics of the PI system. We do this by fitting the same 
statistics for the spectrum of the complete quartz cube to Eq. ( p^ ) with / = 1 
and obtain q = 0.082. We then compare Eq. (|l9|) with the spectral rigidities 
of the other spectra considering / as a free parameter. The best-fit values of / 
for spectra corresponding to different radii of the removed octants are given in 
Table 1. We use these values to calculate the NNS distributions. The results 
of the calculation are compared with the experimental data in Fig. 4. We 
see that the proposed model presents a satisfactory description of the whole 
transition by varying a single parameter. The only disagreement is between 
the calculated and measured values of P{s) at small values of s as we have 
already expected. The symmetry breaking decreases the probability of finding 
degenerate levels sharply leading to the observed dip at small s in the spacing 
distributions |^ . This dip is followed by an overshoot to restore normalisation. 
The width of this dip, which is about 1/10 of the mean level spacing, provides an 
estimate for the minimum level-splitting below which two levels can be regarded 
as approximately degenerate. 
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4 Summary and conclusions 



Conventional models for symmetry breaking in chaotic systems assume that the 
non-diagonal matrix elements of the symmetry breaking perturbation are statis- 
tically equivalent. They predict a monotonous transition from a mixed statistics 
to the GOE as one switches on and increases the symmetry-breaking pertur- 
bation. This is not the behaviour of the acoustic resonance spectra of quartz 
measured by EUegaard et al We propose another scenario for symmetry 
breaking in which the levels are approximately separated into two groups, one 
with and one without symmetry. This happens when the expansion of any of the 
wavefunctions in terms of the eigenfunctions of the symmetry operator is either 
dominated by a single term or composed of statistically equivalent contributions 
from all its component. Considering a system in which a symmetry that has 
two representations is conserved, we represent the spectrum as a superposition 
of two independent sequences, one for each symmetry representation. Now, as 
we switch on the symmetry breaking perturbation, we remove a number of levels 
from these two level sequences into a third sequence. The latter corresponds to 
the strongly perturbed states in which the symmetry has virtually disappeared. 
Because we have increased the number of spectral partitions, the spectrum will 
look as if it became more regular by the small violation of the symmetry. The 
apparent regularity will continue until the three sequences become equally pop- 
ulated. The further increase of symmetry breaking will reverse the evolution 
of the shapes of the spectral statistics toward the GOE shapes. Namely this is 
the behaviour found by inspection of the acoustic resonance spectra of quartz 
measured by EUegaard et al We see, for instance, from Fig. 4 that it is pos- 
sible to describe the evolution of the NNS distribution and the spectral rigidity 
during the transition from a fully conserved symmetry to a fully violated one 
by varying a single parameter, measuring the fractional density of states which 
are practically unaffected by symmetry violation. 
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Table 1. Best fit values for the fraction / of states not affected as a result of 
violating the point-group symmetry of a quartz rectangular block by removing 
spherical octants of various radii r. The values of / are evaluated by fitting the 
spectral rigidities A3 {L) of Fig. 4. to Eq. (pj[) 



r 





0.5 mm 


0.8 mm 


1.1 mm 


1.4 mm 


1.7 mm 


huge dcfocusing 


/ 


1 


0.49 


0.37 


0.28 


0.23 


0.14 






Figure Caption 



Figure 1. The monotonous dependence of (a) the parameter f3 in the Brody 
distribution (Eq.(||)) for the NNS distributions /'b(s) and (b) the parameter qsv 
in the Seligman and Verbaarschot expression (Eq.(pT|)) for the spectral rigidity 
A3 (L) that interpolate between the Poisson and GOE statistics, obtained by a 
numerical experiment on a chaotic system undergoing a gradual breaking of a 
symmetry according to the Guhr-Weidenmiiller model ||^. 

Figure 2. The spectral rigidity A3 (L) for a chaotic system undergoing 
a gradual breaking of a symmetry according to the scenario proposed in the 
text. The results of the calculation when the fraction of states remaining nearly 
symmetry invariant / = 1, 0.9, 0.667, 0.333, 0.1 and are given by the curves 
labelled by A, B, C, D, E and F, respectively, while the curve P corresponds to 
a regular spectrum. 

Figure 3. The spectral rigidity A3 (L) measured in for the quartz block 
when the symmetry is present (open circles) , and when it is violated by removing 
a tiny octant of radius r = 0.5 mm (closed circles). The straight line labelled P 
is for an integrable system. 

Figure 4. The NNS distributions P{s) and spectral rigidities A3 (L) for 
the different radii of the octant removed from a quartz block: (a) r = 0, the 
flip symmetry is fully conserved, (b) r = 0.5 mm, (c) r — 0.8 mm, (d) r = 1.1 
mm, (e) r = 1.4 mm, (f) r — 1.7 mm, (x) the block with the huge defocusing 
structure. The data are from Ref. H, while the curves are the theoretical 
curves using Eqs. (18) and (|l9|). Values of the parameter / that measures the 



fractional densities of states that keep the symmetry are given in Table 1. 
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